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Описан подход к решению задач исследования многоэкстремальных зависимостей для поиска нескольких 
экстремумов одновременно, таких, как определение нескольких равнозначных глобальных экстремумов 
и выделение, наряду с глобальным, наиболее значимых локальных экстремумов. Показана эффективность 
применения для решения таких задач модифицированного метода роящихся частиц. 
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Введение. Множество функций, которые соответствуют пространствам поиска реальных практи- 
ческих задач, являются сложными, т.е. имеют определенные особенности, которые усложняют 
поиск экстремумов. К таким особенностям относят: множество равнозначных экстремумов; нали- 
чие большого медленно убывающего плато, что затрудняет нахождение глобального экстремума; 
множество локальных экстремумов и т.д. Для проверки эффективности поисковых методов были 
построены тестовые функции, которые включают одну из перечисленных особенностей, и при 
этом имеются подробные описания их графиков и известные значения экстремумов. К ним, на- 
пример, относят: функцию Химмельблау, функцию Розенброка, функцию Растригина, функцию 
Нейдорфа и т.д. 

Постановка задачи. Большинство современных методов поиска экстремумов ориентированы на 
одноэкстремальные функции. Поэтому при поиске минимумов и максимумов многоэкстремальных 
функций приходится создавать дополнительные алгоритмы, ориентированные на определенную 
функцию. Таким образом, для исследования многоэкстремальных функций необходим обобщен- 
ный метод, позволяющий находить либо все экстремумы, либо группы экстремумов, удовлетво- 
ряющих определенным условиям. 

Среди эвристических алгоритмов был выбран метод, наиболее подходящий для решения 
координатных экстремальных задач, а именно, канонический метод роящихся частиц (КМРЧ). Ис- 
следователи, основываясь на законах механики движения, предложили его модификацию для 
приближения предлагаемой математической модели к описанию поведения роя в природе (учи- 
тывая биологические и социальные признаки) и алгоритм, позволяющий исследовать многоэкс- 
тремальные поверхности [1]. 

Метод делящихся роев (МДР). В КМРЧ частица в каждый момент времени, т.е. на всех итера- 


циях расчета 1, характеризуется вектором координат Хр определяющими ее положение, векто- 


рами скорости Г;, и ускорения [2]. Далее в каждой точке х‚, вычисляется значение целевой 


1 


функции (критерия оптимизации) О.) (где г- номер частицы). Исходя из полученного значения 


(в) по заданным правилам, частица меняет направление и значение ускорения, что влечет изме- 
нение скорости и положения. Закон изменения скорости частицы выражается в виде уравнения: 


Ь 
Г, Е о; + 019, В] ® (Хх: № х; 1) , (1) 
где х” - вектор координат частицы с наилучшими значениями целевой функции О(), заданной 


в виде функциональной зависимости (4х); И [0,В] - вектор псевдослучайных чисел, рав- 

номерно распределенных в интервале [0, В]; а — свободный параметр, определяющий инер- 

ционное свойство частицы; У,, — скорость 1-й частицы в 1-й момент условного времени. 
Закон изменения положения частицы выражается в виде уравнения: 


Хар = МЕТ Ин (2) 
С целью применения КМРЧ к многоэкстремальным функциям выполняется его модифика- 
ция в МДР. Первая модификация связана с введением понятия «антисоциальной» частицы. Такие 


частицы двигаются в противоположную сторону от центра притяжения, что позволяет обнару- 
жить другие локальные и глобальные экстремумы, иными словами, это частица глобального поис- 
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ка, повышающая эффективность обнаружения экстремумов и снижающая риск зацикливания на 
локальном экстремуме. 

Вторая модификация связана с организацией параллельного роения нескольких роев. При 
обнаружении «антисоциальной» частицей области предположительного экстремума она становит- 
ся центром нового роя, т.е. запускается параллельный процесс роения. Обычно новый рой имеет 
такие же характеристики, как и предыдущий. 


В МДР рой характеризуется двумя радиусами: исходный построенного роя В? и мини- 


мальный, т.е. создается не центр, а область притяжения роя на некотором расстоянии и" от 





точки с координатами частицы. 
Введенные модификации нашли отображение в математической модели КМРЧ. Уравнение 
закона изменения положения частицы (2): 
| ей |. № 
а в ты те о 1 : К 1, К) / (3), 
где К, —- количество роев, образовавшихся к {-у этапу; Л — темп поиска, введенный в рассмотре- 
ние для облегчения настройки процесса поиска; А — номер роя. 
Уравнение изменения скорости (1) принимает вид 
| | |. = 
и =. Е И (4), 


1.1 1,11 








‚Л > 
где А — общее ускорение частицы, которое определяется действием различных сил (сила притя- 
жения и торможения): 
А ВАР + А, (5) 
где ц — коэффициент «социальности» поведения частицы, и = +1 -— для «социальных» частиц, 
т.е. частиц локального поиска, и и = -—1 — для «антисоциальных» частиц, т.е. частиц гло- 
бального поиска; р — коэффициент силы притяжения; {" — коэффициент силы торможения. 
В формуле (5) каждая составляющая ускорения рассчитывается с учетом флуктуирующего 
параметра на основе случайной функции с симметричным распределением относительно номина- 
ла его настройки 


г(^,=)=^| 1+2=(та()-0,5) |, (6) 


где /), - номинальное значение флуктуирующего параметра; = -— отклонение от номинального 
значения. 
Таким образом, результирующее уравнение закона изменения скорости для частицы ло- 
кального поиска (ЧЛП) и частицы глобального поиска (ЧГП) имеет вид 


ко И Кехи ре ой р 
а > Ри = | К )* (^. ‚Е. ) а и (^, р = р й (%.›5, ) Ч (7) 


где г.(^..,=,.) - обеспечивает неоднозначность положения центра притяжения; х,(^„,=»„) — обеспе- 


чивает неоднозначность величины ускорения; г„(.,=) — обеспечивает неоднозначность ве- 


К,ехи 


личины коэффициента трения; х;’’ -— /-я координата наилучшего по критерию О поло- 


ке 


жения точек д -го роя за всю историю его движения от 0 до 1; В. р 


- проекция отрезка 


К, ехи 


вектора х’,’ -х „‚ ДЛИНОЙ и" на 7-ю ось параметрического пространства поиска. 





В результате роения исходного роя по законам изменения (3) и (7) ЧЛП будет найден 
один из экстремумов. В процессе роения заданный процент ЧЛП-частиц на определенное время 
становятся ЧГП, что позволяет осуществить поиск областей предположительного нахождения 
других экстремумов и построить новые рои, которые начнут свое параллельное роение независи- 
мо от стартового. В результате параллельно с исходным может роиться любое количество роев, 
каждый из которых может либо найти новый экстремум, либо сместиться в зону роения другого. 
Это приводит к объединению обоих роев в один. 
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Таким образом, МДР позволяет выполнить эффективный поиск множества экстремумов 

сложных функций. Результативность МДР иллюстрируется примерами его применения для нахож- 
дения экстремумов тестовых функций. 
Примеры применения МДР для поиска экстремумов сложных функций. Особенность всех 
эвристических алгоритмов заключается в невозможности использования математического иссле- 
дования их свойств и выработки эффективных настроек. В связи с этим исследование таких алго- 
ритмов осуществляется методами имитационного моделирования на статистически значимых вы- 
борках вариантов объектов оптимизации. Таким образом, исследование МДР происходит в про- 
цессе его применения для решения задачи поиска экстремумов усложненной одномерной функ- 
ции Нейдорфа и двумерной функции Химмельблау на заданном диапазоне. 

При исследовании функций для одного опыта строится 50 случайно сформированных роев 


с заданным начальным радиусом К, и различными центрами. Процесс роения для каждого роя 
происходит за заданное число итераций №, № =200. 


Метод роящихся частиц с механизмом деления роя, как и большинство эвристических ал- 
горитмов, основывается на множестве параметров настройки, т.е. результативность метода зави- 
сит от выбранных значений этих параметров. При этом значения некоторых из этих параметров 
зафиксированы с расчетом на правильность выбранных значений, с учетом оценок, полученных 
из предыдущих опытов МРЧ. К таким параметрам относятся: деморазмер роя (т), исходный ради- 


ус роя (К, ), флуктуирующие составляющие ускорения (^.,, =), флуктуирующие составляющие 


р 
трения (^.,, =,), флуктуирующие составляющие центра роя (^.,, =.). Кроме этого интуитивно 
понятно, что данные параметры не являются значимыми для МДР. В свою очередь, значения па- 
раметров — конечный радиус А,, процент ЧГП (№...) и число итераций, на которые частица из 
группы ЧЛП переходит в группу ЧГП (Н ), варьируются, так как они напрямую связаны с объек- 
том исследования и являются ключевыми в МДР. Поэтому для оценки влияния ^,, №, иН на 
результирующее решение по трем критериям и совместного влияния всех параметров проводится 
полнофакторный эксперимент (ПФЭ). Исследование по критериям оценки проводится по резуль- 
татам роения 50 роев (№ =50), состоит из 5 независимых опытов для каждого из 8 различных 
условий. В итоге необходимо провести 2000 опытов, что позволит оценить устойчивость резуль- 
татов и влияние случайных исходных данных на эффективность оптимизации. 


При проведении ПФЭ оценка свойств МДР, по опыту проведенных исследований и сути 
решаемой проблемы, проводится по четырем важным и показателям: 


-Е (1; }), эффективность нахождения всех точек экстремума исследуемой области. Дан- 
ный критерий характеризует метод применительно к многоэкстремальным задачам. Он формиру- 
ется как относительное среднее количество экстремумов, найденных МДР, заранее известному 
количеству экстремумов исследуемой функции; 

_5 (1; }), ресурс, который определяет минимальное количество итераций, необходимых 
для нахождения локального экстремума. Для данного исследования он вычисляется как среднее 
число итераций, за которое центр роя достигает области нахождения одного экстремума, относи- 


тельно общего числа итераций локального поиска, № = 200; 

-А(}а}), точность оценки величины экстремума, т.е. средняя оценка близости найденно- 
го значения оптимума к заранее известному при тестировании значению экстремума; 

- 5е (<. !), глобальный ресурс, определяет количество итераций, за которое будут найде- 


ны все экстремумы. Он вычисляется как среднее количество итераций, за которое центры множе- 
ства роев достигли области нахождения всех экстремумов к числу найденных экстремумов. 

При этом эффективность метода является ключевым критерием, а три остальных - вспо- 
могательными и могут накладывать дополнительные ограничения на значения параметров, если 
они не противоречат результатам оценки по критерию А. 
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Пример 1. Исследование многоэкстремальной одномерной функции. 
Для исследования выбрана следующая функция Нейдорфа, график которой изображен на 
рис.1: 


_ 1-х)" +15зш(Ых +2) 


(8) 
\ х° + 0,5 


Исследование функции проводится на интервале хе |[-20,20], где функция имеет 6 мини- 
мумов: (-13,72; 15,22), (-8,38; 9,03), (-2,83; 0,42), (2,04; —5,7); (7,64; 4,12) и (13,01; 10,42). 


У(х) 


У(Х) 40 


30 


20 


10 





Рис.1. График функции Нейдорфа 


Ставится задача нахождения экстремумов рассматриваемой функцией МДР. Для решения 
данной задачи были выбраны следующие параметры эксперимента: т = 50 частиц; К, =2 


(5 % от размера пространства поиска); ^.„ =0,06, =, =1; Л, =1, в, =0,7; |, =1, 0,6, 0,2; 


№, =10%, 20%, 30%; Н =5,15,20. 


Анализ результатов проведенного ПФЭ показывает, что для оптимальной эффективности 
МДР, т.е. нахождения всех экстремумов, необходимо установить максимальное значение пара- 
метров ^, иН. В свою очередь, оптимальные значения по другим критериям оценки получаются 





при минимальном значении параметра №... Таким образом, в качестве эффективных значений 
факторов, при которых достигаются лучшие значения по всем критериям, выбираются: К, =1, 
№, =10 иН = 20. 

Для уточнения полученных значений параметров настройки проводится дополнительный 


ПФЭ с более узкими интервалами варьирования факторов. В качестве центра плана выбирается 
точка эффективных значений из предыдущего ПФЭ. 

В результате происходит уточнение значений параметров настройки: А, =1, №... =10 
и Н = 25. При данном наборе факторов критерии оценки имеют следующие значения: Е=1 (най- 
дены все 6 экстремумов), $ = 0,043, А = 2,77Е-07, $4 = 0,0304. 

Процесс выполнения МДР проиллюстрирован на рис.2. По оси Х откладываются итерации 
роения, а по оси У — значения функции в центре роя на каждой итерации. Сплошными линиями 
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обозначены траектории движения центров роев, которые в итоге нашли экстремумы, пунктирны- 
ми линиями, в свою очередь, обозначены траектории роев, которые перешли в область локализа- 
ции уже найденных экстремумов. Кроме этого представлена иерархия роев, т.е. показано, сума- 


сшедшие частицы какого роя создали следующие рои. 


Е 


Для более наглядного представления траекторий движения роев, нашедших экстремумы, 





{(-28,168)=30,54/ 


я 
‚ {-13,723)=15,224 


#(13,032)=10,418 `.  113,008)=10,418 


=(-3,376)=9,027 = 
ь {(-8.382)=9,027 
Е(7.642)=4,124 


{(-2,752)=0,436 32)=0,416 


1-2; 
о >> 
со 09 
— — 


#(2,039)=-5,698 


Рис.2. Графическое отображение процесса МДР 


они строятся в логарифмических координатах по значению функции (рис.3). 


значение функции 


Г мм 


я 





1, = 
т, 
#2 


МН ГИ = 
Мн ГИ = 
М Г 
а г 
а гг 
М == 


ЧС ЕВЕНЕНЕНЕЬЕЬ 
ЕЕ 


3.2515. 
С, 1525 
3,1553 


ОРЕЛ ЕЕ 


— ЗЕ 
—45, ЗЕ 
И 
= 
ЗЕ 





—, ЗЕЕ 


АС, УЕРЕ 
О, те 
Я, УЧЕНЕЕ 
и 


Пи =. НТ = рееааални в 
Рис. 3. Траектория движения роев в логарифмических координатах 
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Таким образом, показана эффективность применения МДР для нахождения экстремумов 
одномерной функции при определенных значениях параметров настройки. 
Пример 2. Исследование двумерной функции Химмельблау (ФХ). Следующая по сложно- 
сти тестовая многоэкстремальная функция Химмельблау. Она имеет четыре глобальных миниму- 
ма, значения в которых равны нулю: (-2,8051; 3,1313), (-3,7793; -3,2832), (3,2), (3,5844; -1,8481). 
На рис.4 представлен график поверхности и линии уровня функции Химмельблау с отмеченными 
экстремумами. 





(перевернутый вид) 


Рис.4. График поверхности и линии уровня функции Химмельблау 


Для решения данной задачи выбираются следующие параметры эксперимента: 

— фиксированные параметры: т=800 частиц, К, =6, Л,„=0,5, =, =0,05, &, =0,071, 
А, =0,71, А,, =1, &, = 0,08; 

— варьируемые параметры: Л... = 10 %, 20 %, 30 %, 


С 





К.|=1,5;1;0,5, Н =5,15,20. 

В результате проведения серии ПФЭ получается вектор эффективных значений факторов: 
К, =1,5, №... =30 и Н = 20. При этом получены следующие оценки критериев: Е = 1 (найдены 
все 4 экстремума), $ = 0,025, А = 3,944, Е = 31, $4 = 0,0071. 


значение функции 


Е резщеры! 





Рис.5. Процесс МДР для ФХ 


Отображение процесса МДР представлено на рис.5. В связи с тем, что все экстремумы имеют ну- 
левое значение функционала, то все центры роев находят одно и то же значение функции, но 
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при различных значениях аргументов. В данном случае рой 1, 2, 4 и 5 нашли глобальные мини- 
мумы, а рой 3 переместился в область локализации экстремума, который был найден 
роем 2. Для большей наглядности процесс деления роев представляется в логарифмических ко- 
ординатах по оси «Значение функции» (рис.6). 


= 





Рис.6. Процесс МДР для ФХ (в логарифмических координатах) 


Заключение. Показано, что метод деления роя является эффективным инструментов для поиска 
множества экстремумов сложных функций, т.е. для решения многоэкстремальных задач. Эффек- 
тивность МДР зависит от выбранных значений параметров настройки, которые могут отличаться 
при решении различных задач. 
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